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Abstract. 

We study the problem of how a "living" system complex in structure can respond perfectly to local 
changes in the environment. Such a system is assumed to consist of a distributed "living" medium and 
a hierarchical "supplying" network that provides this medium with "nutritious" products. Because of 
the hierarchical organization each element of the supplying network has to behave in a self-consistent 
way for the system can adapt to changes in the environment. 

We propose a cooperative mechanism of self-regulation by which the system as a whole can react 
perfectly. This mechanism is based on an individual response of each element to the corresponding 
small piece of the information on the state of the "living" medium. The conservation of flux through 
the supplying network gives rise to a certain processing of information and the self-consistent behavior 
of the elements, leading to the perfect self-regulation. The corresponding equations governing the 
"living" medium state are obtained. 

Key words, active hierarchical systems, cooperative self-regulation, information self-processing, 
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1. Introduction. "Living" systems and the self-regulation problem. 

The present paper is devoted to one of the fundamental problems of how complex 
"living" systems widely met in nature can adapt to changes in the environment. By 
the term "living" system we mean one that comes into being, provides for itself, and 
develops pursuing its own goals. This class comprises a great variety of biological 
and ecological systems. Besides, economic systems can be also placed into such a 
class because their occurrence, growth, and development is due to self-organization 
processes. To make the subject of our analysis more clear let us, first, specify the 
main properties of the systems we deal with. 

1.1. Living systems. The elements of such a system, first, should be perma- 
nently supplied with external "nutritious" products for their life activities. Second, 
in order for its elements to subsist special conditions are required. The latter arc 
implemented through these elements making up a certain medium (which will be re- 
ferred below as to a "living" medium) whose state is controlled by keeping its basic 
parameters inside certain vital intervals. 

As a rule both of these requirements are fulfilled by a supplying network which 
provides the elements with nutrients as well as controls the living medium state. Since 
the external products usually penetrate into the system through a common entryway 
and, then, should be delivered to a great number of elements the flow of these products 
has to branch many times until it reaches these elements. So the supplying network 
is to be organized hierarchically and to involve many levels. 

Let us, now, present typical examples of such systems. First, this is living tissue, 
where blood flowing through a vascular network supplies cells with oxygen, nutritious 
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products, etc. At the same time blood withdraws carbon dioxide and other products 
resulting from the cell life activities, keeping their concentrations inside the vital 
intervals (see, e.g. jij). In this way blood flow controls also the tissue temperature 
We note that the temperature and the concentration of carbon dioxide (or a 
similar substance) are major parameters characterizing the living tissue state 0, 
because their values are directly determined by the life activities of cells. In order to 
supply with blood each small group of cells the vascular network or, more precisely, 
its arterial and venous beds are approximately of the tree form and involve many 
vessels of various lengths and diameters. The regional self-regulation processes keep 
the carbon dioxide concentration and the temperature inside certain vital intervals 
only in which living tissue can function normally (see, e.g. |l], 

A similar example is a respiratory system where oxygen going through a hierar- 
chical network of bronchial tubes reaches small capillaries. 

Second, large firms are a clear example of economic hierarchical systems. Man- 
agers of all functions and of all levels make up a management network || . Roughly 
speaking, the management network controls both the money flow toward the orga- 
nization low level comprising workers and also the flow of products in the opposite 
direction. In performing technological processes the wages paid to workers actually 
transform into the firm products. 

The existence of tremendous amount of goods in the market, in contrast to a 
relatively small number of raw materials shows that there also must be large hierar- 
chical systems supplying the consumers with goods. The goods flow after reaching the 
consumers transforms into money flow in the opposite direction [0, 0]. In particular, 
self-organization phenomena can give rise to trade networks of the tree form Q . 

Concerning ecological systems we would like to note that they are also complex in 
structure and can involve a larger number of "predator-prey" levels . The dynamics 
of ecosystems is governed by biomass and energy flow on trophic networks and under 
certain conditions leads to formation of the tree like trophic networks ||] . 

Keeping in mind the aforementioned examples let us represent a living system 
as one consisting of a living medium and a supplying network which involves sup- 
plying and draining beds of the tree form (Fig. 1.1). (In principle, the two beds can 
coincide with each other in space.) A transport agent (e.g., blood in living tissue) 
"flowing" through this network delivers "nutrients" to the living medium and with- 
draws "life activity" products. The interaction between the transport agent and the 
living medium (leading to the product exchange) takes place when the transport agent 
"flows" through the last level branches of the supplying network (in Fig. 1.1 through 
the branches of level N) . 

The concentration of the life activity products will be regarded as the parameter 
characterizing the living medium state. This is due to the fact that the higher are 
life activities in intensity, the greater is the concentration of the resulting products. 
So an increase in this concentration informs the system that the living medium needs 
a greater amount of "nutrients" . In addition, the products resulting from the life 
activities and remaining in the living medium themselves can depress its functions. 
Therefore the system has to prevent the concentration of these products exceeding a 
certain critical value. 

The motion of transport agent is accompanied by energy dissipation. So a certain 
external force should be applied to the system that affects the overall flow of trans- 
port agent (for example, in living tissue this is the blood pressure Jj], ||). Further 
distribution of the transport agent flow over the supplying network is governed by 
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Fig. 1.1. "Living" system structure. (For simplicity is shown the binary network.) 

the individual properties of its branches. More precisely, the transport agent flow 
distribution is controlled by both a certain potential (in living tissue it is the blood 
pressure) and by the "resistances" of the branches to the transport agent flow. In 
order to specify the physical regularities of the transport agent flow distribution we 
may choose various approaches. One of them is Prigogine's principle of minimum 
entropy production stated in nonequilibrium thermodynamics (see, e.g., Q]). An- 
other is related to Swenson's principle of maximal dissipation. In our analysis the two 
approaches lead to the same result: the transport agent flows meet an extremum of a 
certain functional. So, to be defined we will adopt the minimum entropy production 
principle. 

1.2. Self-regulation problem. For different elements of the living medium to 
fulfill their individual functions independently of one another the supplying network 
should, at least in the ideal case, have a capacity for controlling the state of the 
living medium at each its point. This requirement is reduced to the local control 
of the "perfusion" rate of the transport agent through the living medium. So the 
independence of the perfusion rate at one point from life activities of elements located 
at other points is a desirable property which will be referred below as to the perfect 
self-regulation. In fact, if we ignore diffusion in the living medium solely the perfusion 
rate of transport agent at a given point will govern functioning of elements at this 
point: supply them with "nutrients" and withdraw the life activity products. So 
the perfusion rate should take such a value that provides the optimal conditions 
for functioning of the elements, depending on their life activities. Therefore if the 
perfusion rate at one point reacted substantially to the life activities at other points 
then the living medium elements would interfere with one another and the living 
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system could lose its capacity for adapting. The existence of a great variety of living 
systems in nature enables us to think that such a nonlocal interaction is suppressed 
or, at least, depressed to a certain degree. However, in trying to describe how this self- 
regulation can be implemented, i.e. its particular mechanism we meet the following 
fundamental problem. 

The perfusion rate of the transport agent at a given point r is certain to be deter- 
mined by its flow through the corresponding last level branch zjv(r) of the supplying 



network (Fig. 1.1). However, this flow in turn depends on the transport agent flow 
through the branch ijv-i(r) of the previous level connected with the branch ?at(i") 
and so on up to the stem io (the branch of zeroth level). In this way we find a path 
P(r) = {io, «i(r), . . . , ijv-i(r)> ijv(r)} on the supplying network through which the 
transport agent flow reaches a small neighborhood of the point r. Let us consider the 
transport agent flow through one of these branches, for example, a branch i n (r) of 



level n < N (in Fig. 1.1 it is the branch ?Ar_2(r)). This branch supplies with "nu- 
trients" not only elements in the vicinity of the point r but also all of the elements 
located inside a certain domain 9Jl„(r) (inside the domain 9Jtjv-i(r) in Fig. 1.1) and 
controls the living medium state in this domain as a whole. So the required flow of 
the transport agent through the branch i n (r) is specified by the life activities of all 
the elements belonging to the domain 9Jl n (r) and a change in the life activities of one 
of them inevitably will cause this flow to vary. 

Since the external force applied to the system controls only the overall flow the 
specific distribution of the transport agent flow over the supplying network is governed 
by the "resistances" of all the branches. Under such conditions, in general, a change in 
the transport agent flow through the branch i n (r) (for example, because of variations 
in its "resistance" ) causes the transport agent flow to alter in all the branches of higher 
levels that are connected with i n (r), in particular, in a last level branch i./v(r') leading 
to a point r' 6 9Jt„(r) and r' ^ r. So, in principle, a change in the life activities of 
elements near the point r will lead via the branch i n (r) to variations in the perfusion 
rate of transport agent at all the points belonging to the domain 9Jt„(r). Since the 
domains {9Jt„(r)} increase in size and tend to the total space of the living medium 
as we pass on the supplying bed from the last level to its stem this effect is nonlocal 
substantially. 

Therefore special conditions are required for the self-regulation to be perfect be- 
cause in the general case the perfusion rate of the transport agent is determined by 
the life activities of all the elements. 

In other words for the perfect self-regulation to occur the "resistances" of all the 
branches should vary in a self-consistent way. In this place we actually meet the 
fundament problem mentioned above: what governs such variations of the branch 
"resistances" ? For a large living system it is unlikely that any one of its elements can 
possess the whole information required of governing the perfect response to changes 
in the environment and, so, in the life activities of living medium. This is due to 
the fact that such a control requires processing of a great amount of information 
characterizing the living medium state on all spatial scales (i.e. at all the levels of 
supplying network). In particular, in living organisms local alterations in homeostasis 
seem to be controlled by regional mechanisms of self-regulation rather than the central 
nervous system Jl|. Besides, none of these elements can individually control transport 
agent flow through the supplying network because of the mass conservation at its 
branching points. 

The aforementioned allows us to assume that the perfect self-regulation can be 
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implemented through cooperative mechanisms. By this scenario each branch of the 
supplying network receives the corresponding small piece of the information on the 
living medium state and their reaction to this information gives rise to the desired 
redistribution of the transport agent flow over the supplying network. In other words, 
the branch "resistances" vary in such a self-consistent way that enables the supplying 
network to provide, for example, an additional amount of "nutrients" for the living 
medium elements which have a need for this and not to disturb transport agent flow 
at other points of the living medium. 

However, for the cooperative mechanism of perfect self-regulation to come into 
being, first, a certain self-processing of the information is required which enables the 
branches to react adequately, depending on their place in the supplying network. 
Second, the physical properties of the supplying network should give rise to the co- 
operative effect of variations in the branch "resistances" on the transport agent flow 
redistribution leading to the desired results. 

The purpose of the present paper is to show that these requirements can be ful- 
filled and so a cooperative mechanism of perfect self-regulation in living systems can 
exist. We will demonstrate that there is a required self-processing of information and 
the physics of transport agent motion ensures its proper redistribution over the sup- 
plying network. Beforehand we can say that it is the mass conservation at the nodes 
of supplying network and its hierarchical organization those bring into being the in- 
formation self-processing. This information self-processing is implemented through 
measuring the concentration of the life activity products inside each branch of the 
draining bed. Due to the interaction with the living medium the transport agent 
saturates with the life activity products when it flows through the last level branches. 
Then the transport agent moves through the draining bed from higher hierarchy lev- 
els to lower ones and, so, at the nodes its smaller streams unite into larger streams. 
Therefore, due to the mass conservation the concentration Qi of life activity products 
inside, for example, a branch i is actually equal to the concentration of life activity 
products inside the living medium averaged over the domain 9Jti which is drained as 
a whole by this branch. So the value of Qi aggregates the information on the living 
medium state in the domain dJli as a whole and, thus, can play the role of the infor- 
mation piece needed for this branch and the corresponding branch of the supplying 
bed to react properly to variations in the life activities of the living medium. The cor- 
respondence between the branches of the supplying and draining beds is illustratied 



in Fig. pL.l| by the pair i N -2{r), i* N ^ 2 {r). 

It should be noted that developing this model we actually have kept in mind 
one of the possible mechanisms of self-regulation in living tissues. Life activities of 
cellular tissue gives rise to variations in the carbon dioxide concentration in it and, so, 
in venous blood. Receptors embedded into the vein walls through a regional nervous 
system govern the expansion or contraction of the corresponding arteries (see, e.g., 
il). 

Below, at first, we will analyze in detail the model for a living system with a 
regular supplying network. This network is organized in such a manner that the 
perfusion rate of transport agent be the same at every point of the living medium, 
all other factors being equal. In particular, we will develop a technique that not only 
enables us to study the perfect response of living system to changes in the living 
medium state but also can be used in investigations of its more complex behavior. 
Then we will show that the perfect self-regulation can also occur in living systems 
with nonregular supplying networks of the general form. 
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2. Mathematical model for living system with regular supplying net- 
work. Let us consider a system consisting of a living medium 971 and a hierarchical 



supplying network 91 similar to that shown in Fig. 1.1. The network 91 involves sup- 
plying (91 s ) and draining (91 d ) beds of the tree form (illustrated in Fig. 1.1 by the 
left and right-hand side parts) and an external element £ C xt joining the tree stems. 
Transport agent flowing through the former bed supplies the living medium 971 with 
"nutrients" . At the same time transport agent withdraws life activity products from 
the living medium through the draining bed. For simplicity we assume that the living 
medium domain 9Jt is a d-dimensional cube of edge Iq and the geometry of both the 
beds is the same. 

The supplying and draining beds can be represented as the collections of groups 
of branches {i} belonging to one level: 91 s = {L^I^Lq, 9T d = {L^}^Lo> where 
h s n {d) = and i e level n). Below we will use the symbol L„ instead 

of or L d where it does not lead to misunderstanding or gives the same results 
due to the mirror symmetry of the supplying and draining beds. In these terms the 
embedding of the supplying network 91 (as the collection of branches {?}) into the 
living medium matches the partition of the domain 9Jt into the following collection of 
cubes {{OT} ieL „}£U. 

For each level n the group M„ = {dJl 1 }^^ involves 2 nd equal disjoint cubes of 
edge l n — Iq2~ u that together compose the whole living medium domain: 



(2.1) 9Jt=(j9Jr and 07T f] TV = if i ^ j, i,j G L„. 

iGL„ 

For any two levels n and m < n each domain 971* £ M„ belongs to one of the domains 
of the group M m and, thus, each domain 9JP G M m can be represented as a certain 
union of the domains of group M n 



(2.2) 97P = |J an*. 

For zeroth level 9JT° = 971. Below the cubes of the group M„ = {97l J } ieL „ will be 
also called fundamental domains of level n and those of the last level N will be also 
referred to as elementary domains. 

The supplying and draining beds are embedded into the living medium in such 
a way that every branch i of a given level n supplies (or drains) a certain domain 
97T G M n as a whole. Each elementary domain is bound up with one of the last level 
branches. The last level number N is assumed to be much larger than unity: N ^> 1, 
and the length In may be regarded as an infinitely small spatial scale. Transport 
agent flow is directed from lower to higher levels on the supplying bed and in the 
opposite direction on the draining bed. 

Since the geometry of both the beds is assumed to be the same, we specify it for 



the supplying bed only. The stem of this bed (Fig. 2.1) goes into the cube 971 through 
one of its corners and reaches the cube center Oq, where it splits into g — 2 d branches 
of the first level. Each branch of the first level reaches a center 0\ of one of the g 
fundamental domains of the first level. At the centers {Oi} each of the first level 
branches in turn splits into g second level branches. Then this process is continued in 
a similar way up to level N. The branches of the last level N are directly connected 
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Fig. 2.1. Fragments of the supplying network architectonics for the two- and three-dimensional 
living medium 

with living medium 9Jt. It should be noted that in this model the length £ n of branches 
belonging to level n is estimated by the expression £ n <~ \fdl n = \[d2~ n and smaller 
is a branch, to a higher level belongs it. 

Transport agent does not interact with the living medium 971 during its motion 
through the branches except for the last level ones. When the transport agent reaches 
one of the last level branches of the supplying bed, for example, branch i n the "nu- 
tritious" products carried with it uniformly spread over the elementary domain TV N 
containing the branch ipj. At the same time transport agent is saturated with the 
life activity products, withdrawing them from the domain TV N . In this way the life 
activity products going with the transport agent through the draining bed leave the 
system. 

The state of the medium 971 is described by a certain field 9(r,t) being the di- 
mensionless concentration of life activity products. Dynamics of the field 0(r,t) is 
governed by the volumetric generation q(r,t) resulting from the life activities, the 
dissipation because of draining the living medium, and the diffusion. In other words 
the field 6(r, t) is considered to evolve according to the equation 



where D is the diffusivity and f](r,t) is the perfusion rate of transport agent. Since a 
branch i N of the last level supplies the elementary domain 9Jt l N OI the living medium 
as a whole the transport agent flow J iN through the branch i N and the perfusion rate 
r](r, t) are related by the expression 



(2.3) 



— = DX7 2 8 + q-8 V , 



(2.4) 




N 



The distribution of the transport agent flow { Ji} over the network 9T obeys the mass 
conservation at its nodes {b} 
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(2.5) 



ieb 



where is the transport agent flow through a branch i going in or out of the node 
b and the sums run over all the branches i € b leading from or to this node for 
the supplying and draining beds, respectively. Provided the field r)(r,t) is given 
expressions (2.4) and (2.5) completely specify the quantities {Ji}. 

The last term in equation (2.3) implies that the transport agent when going 
through the branch is saturated with the life activity products up to their concen- 
tration in the elementary domain OKJv- Then these products moves with the transport 
agent through the draining bed without exchange with the living medium. Therefore, 
we also ascribe to the transport agent going through the draining bed the collection 
of variables {9i} measuring the concentration of the life activity products inside the 
branches {i £ 9T 1 } and regard the value Ji9i as the flow of these products in the 
branch i. For a branch ijv of the last level we write 



(2.6) J lN 9 lN = J 0(r)tj(r)dr. 

and assuming the conservation of the life activity products at the nodes of the draining 
bed we get 



(2.7) 



ieb 



(Ji< 



Equation (2.3) among with expressions (2.4) and (2J3) describe the dynamics of 
the state of the the living medium and its interaction with transport agent. Equations 
(2.5) and (2.7) reflect the general laws of transport phenomena in the supplying and 
draining beds. 

Let us, now, specify the regularities governing the transport agent flow through 
the network 91. First, provided the characteristic parameters of branches are fixed 
the transport agent flow patter {J,} is assumed to meet the minimum condition of 
the total energy dissipation due to transport agent motion: 



(2.8) D{ JA min 

{■/*} 

subject to relationships (^.5j), where the total rate of energy dissipation is given by 
the expression 

(2.9) P{ Ji} = \J2 R ^ ~ Jo£cxt - 

i 

Here Ri is the kinetic coefficient characterizing the energy dissipation during the 
motion of transport agent through the branch i, £ ext is the external force causing the 
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transport agent motion through the network 91 as a whole, Jo is the transport agent 
flow through the tree stems, and the sum runs over all the branches of the network 91. 
It should be noted that previously by the term "branch resistances" we have exactly 
meant the given kinetic coefficients because they do play the role of the true branch 
resistances to transport agent flow as it will be shown below in the next section. 

Second, in this model it is the coefficients {Ri} those characterize the individual 
effect of the branches on the transport agent flow. So the system response to changes 
in the living medium state, namely, to variations in the field 6(r, t) (see Introduction) 
is represented as time variations of the coefficients {Ri}. Since the supplying and 
draining beds have been assumed to be equivalent in architectonics the flow pattern 
on these beds {■/,} as well as the coefficient collection {Ri} are set to be the mirror 
images of each other within reversing the flow direction. The latter allows us to 
confine our description of the system response to the draining bed. 

As mentioned in Introduction the pattern {6^} can be treated as the aggregated 
information on the living medium state on all the spatial scales. In particular, the 
variable 9i corresponding to the branch i characterizes the state of the living medium 
in the fundamental domain 971* as a whole. Therefore we assume that for each branch, 
for example, a branch i of level n time variations in the coefficient Ri are directly 
controlled by the variable 9i assigned to this branch. Under the steady-state conditions 
this means that the coefficient Ri is an explicit function of the variable 9i identical 
for all the branches of one level n, i.e. Ri — R n (6i). The general properties of the 
R n {8) dependence are actually determined by the fact that the transport agent flow 
through the network 91 should grow as the variables {6i} increase. Indeed, let changes 
in the environment cause the life activities of the living medium to grow in intensity. 
The latter leads immediately to an increase in the concentration 9(r, t) of life activity 
products inside the living medium, causing an increase in the variables {9i}- The 
higher is the life activity intensity, the greater is the amount of "nutrients" needed for 
the living medium. Therefore, under such conditions the transport agent flow through 
the network 91 must increase too. In other words, an increase in the variables {6^} 
should give rise to an increase of transport agent flow. The less is the coefficient Ri, 
the greater is the transport agent flow that can go through the branch i, all other 
factors being equal. Thus, the function R n (9) must be decreasing with respect to 
the variable 9. In addition, there should be a certain critical concentration 9 C of the 
life activity products showing the upper boundary of the vital interval for the field 
0(r,t), i.e. the maximum of the allowable concentration of the life activity products 
in the living medium. The system has to prevent the field 6*(r, t) exceeding the value 
9 C as much as it can. Otherwise, the system can lose the capability for adapting. 
In particular, if the field 9(r, t) reaches the boundary of the vital interval at all the 
points of the living medium and, consequently, all the variables {6^} come close to the 
critical value 9 C (as it follows from ( |2.5| )-( |2~7| )) all the branches should exhaust their 
capability for decreasing the coefficients {Ri}. 

Taking the aforesaid into account we represent the R n (9) dependence in terms of 



(2.10) R n (9) = R° n ^ ¥c ), 

where is a constant equal to R n (9) at 9 — and cf>(x) is a certain universal f unc tion 



of x = 9/9 c . The characteristic form of the </>(x) dependence is shown in Fig. |2.2| a. 

For the transport network 91 to be able to react properly to local variations in 
the life activities the transport agent flow should be governed by branches of all the 




Fig. 2.2. Functions describing the branch response. (The dotted lines correspond to the ideal 
case.) 



levels. As will be seen below this condition allows us to represent the dependence of 
the value R® on the level number n in the form 



(2.H) Rl = Rog n P (n), 

where Rq is a certain constant and p(n) is a smooth function of n such that p(0) = 1 
and formally p(n) — > as n — ► oo. 

Concluding the description of the system response we also take into account a 
possible time delay of the branch response to variations in the variables and 
represent the evolution equation for the kinetic coefficient Ri of branch i belonging 
to level n as 



(2.12) 



Ri 



where r n is th e time delay depending solely on the level number n. This form of 
equation ( 2.12 ) enables us to regard the term 8i/0 c as a dimensionless signal gene rated 
by receptors imbedded into the branch i. As follows from ( 2.1C ) and ( 2.12 ) the 
functions f{x) and 4>(x) are related by the expression 



(2.13) [1 - <Kx)]f[4>(x)] = x. 

The behavior of the function f(x) is displayed in Fig. |2.2| b. 

In the context of the stated model the problem of self-regulation is reduced to 
finding an expression relating the perfusion rate ?y(r, t) of the transport agent to the 
field 0(r,t): 



(2.14) ?/ = ${6»} or 77 = ${(9,7?}. 

In general, this expression is of the functional form, i.e. is nonlocal in space and only 
under specific conditions can relate the fields 77 (r, i), 0(r, t) taken at the same point r. 



COOPERATIVE MECHANISM OF SELF-REGULATION 



11 



In what follows we will find this conditions and obtain the corresponding particular 
form of expression ( [2.14 ) 



3. Governing equations for the perfusion rate. General description. By 

virtue of the adopted assumptions the size In of the elementary domains plays the 



role of an infinitesimal length. Therefore expression (2.4) can be rewritten as 



(3.1) 



M = ^ E J ^> 

N iSL N 



where the symbol 0,(r) stands for the characteristic function of the fundamental 
domain 9JL , i.e. 



(3.2) 



9<(r) 



i if re 

if r4m 



and the sum runs over all the branches of the last level. If we find expressions relating 
the pattern {J,} to the fields 8(r,t), rj(r,t) then formula (3.1) will immediately give 
us the specific form of (2.14). So in the present section we solve the equations stated 
above with respect to the variables {Ji} regarding the field 0(r,t) (and may be the 
field r](r,t)) as fixed. 

3. 1. K irchhoff's equations. Let us, first, find the extremal equations for func- 
tional (2^9). Follow ing the Lagrange multiplier met hod we reduce the extremum prob- 
lem for functional (2J)) subject to conditions (2__5) to finding the extremal equations 
of the following functional 



(3.3) 



2 ^ 



RiJi — SextJo 



En 



where {Pb} are the Lagrange multipliers ascribed to the nodes {b} of the network 91, 
the sum J2b runs over all the nodes {b}, and the symbols X^eb ^> Sieb ^ stand for 
the sum over all the branches going in or out of the given node b. (The direction of 
the motion on the network 91 is chosen so to coincide with the direction of transport 
agent flow.) From the conditions dH>L/dJi = 0we obtain the desired extremals 



(3.4) 

and 

(3.5) 



IB— P l 



PI 



pO _ pO _ c 
Mil -'out — °e 



Here the multipliers P* n , P* nt correspond to the entrance and the exit of the given 
branch i and P® n , P° ut are the multipliers ascribed to the entrance and the exit of the 
transport netw ork. It should be noted that equations (3.4) and (3.5) together with 



equations (2.5) actually make up the system of Kirchhoff's equations for the network 
0T. In these equations the variables {Pb} play the role of potentials at the nodes {b} 
causing the transport agent flow through the corresponding branches and the kinetic 
coefficients {Ri} may be regarded as the branch resistances. 
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3.2. Additional potenti al s ou rce s. At the next step we should solve the sys- 
tem of Kirchhoff's equations ( |2.5| ), (3.4), (3.5). However, when the field 0(r,t) is 
nonuniform, the resistances {Ri} of all the branches can differ in magnitude because 
of the network response. In this case solving Kirchhoff's equations directly is trou- 
blesome. In order to avoid this problem we make use of the following trick. 

Let us ascribe to each branch i the quantity e% defined by the expression 



(3.6) 



Si — Ji(R n . — Ri 



where n, is the level number of the given branch. The quantities {si} will be called 
additional potential sources. Formula (3.6) allows us to rewrite equation (3.4) in 
terms of 



(3.7) Ji J R° i =^-P i ut + £i . 

The collection of equations ( |2.5| ) and (3.7) written for every node and every branch 
forms the system of Kirchhoff's equations describing the flow pattern { Ji} on a certain 
network of the same geometry where, however, branches do not respond to variations 
in {9i}. In this case the self-regulation process is effectively implemented through 
the appearance of the additional potential sources {£i}- We will refer to this net- 
work as the homogeneous one. In this way, when the quantities {e{\ have known 
values the analysis of the transport agent flow redistribution over the supplying and 
draining beds with branches sensitive to Qi is reduced to solving the Kirchhoff's equa- 
tions for the corresponding homogeneous network. In the next subsection we will get 
expressions relating the additional potential sources {e^ to the fields 9(r, t), r)(r,t). 
Therefore this trick, in fact, can be useful in solving the original Kirchhoff's equations. 
In the present subsection we obtain the expressions specifying the flow pattern { J{\ 
for fixed values of the quantities {si}. 

Since the flow patterns on the supplying and draining beds are mirror image of 
each other we may confine our consideration to the draining bed only. In this case we 
may set P° ut = at the network exit and P ; l n = £ e xt/2 for all the branches of the last 
level, i £ Ljy ■ For the corresponding draining bed of the homogeneous network the 
solution of Kirchhoff's equations (2.5), (3.7) can be written in the form 



(3.8) 



Ji 



i 



7 A-iio ^ c 



Here Ay is the Green matrix, i.e. the solution of these equations when £ e 







1, 



and io stands for 
results 



and £ji = for all branches except the branch j for which Ej 
the draining bed stem. We note that the possibility of representation (3 
from the linearity of equations ( |2.5] ) and ([T?]) with respect to the transport agent 
flows {Ji}. Appendix |A| specifies the Green matrix Ay at leading order in the small 
parameter \p( n) — p(n + 1)] j pin) (for the definition of p(n ) see expression (2.11)). 

in terms of 



Formula (3^8) enables us to represent relationship 



(3.9) 



W = Td Yl (D A ^ + 2 A «o£cxt)e 4 (r). 

N ieh N j 



This expression will lead us to ( 2.14 ) if we can describe the dynamics of the additional 
potential sources {£i}. 
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3.3. Gover nin g equations for the additional potential sources. Differen- 
tiating relation (3.6) with respect to time t and taking into account equation ( 2.1 2| ) 
we find the following evolution equation for the quantity ef. 



(3.10) 



dt 



dt 



In Ji 



R ni J >.- 



In obtaining (3.1C) we also have assumed that the flow Ji cannot change its direction. 
Due to the conservation laws (2.5) and (2.7) the values Ji and OiJi corresponding to 
the branch i can be directly related to the same variables ascribed to the branches of 
the subtree 9T whose stem is the given branch i. Passing to the last level of this 
subtree we immediately get 



(3.11) 
(3.12) 



Ji=Y.Jj=j drQi(r)r,(T), 



where we have made us of expressions (2.4), fl2.6|) and the identity 



(3.13) 



9 j (r) = 6 i (r) 



which is due to fl2.2|). 

The system of equation ( 3.1 0| ) and expressions ( 3.11 ), ( 3.12|) describes the dy- 
namics of the additional potential sources in such a form that explicitly relates time 
variations of {si} to the fields 0(r,t), r)(r,t) and, thus, among with expression (3.9) 
gives us, in principle, the desired relationship (2.14). 



In particula r, wh en the branches respond without delay, i.e. {r„ = 0}, this system 
and expression ( 2.13 ) gives us the explicit relation 



(3.14) £i(t)=R° drri(T,t)ei(r) 



an 



, / dr7 ? (r,t)^(r,t)9 l (r) 

m 

J dr V (r,t)Qi(r) 
m 



The substitution of ( 3.14 ) into (p.9[) directly specifies the functional dependence (p,14|) 



4. Ideal supplying network. In this section we show that there are special 
conditions under which the supplying network can function perfectly, namely, local 
variations in the concentration of life activity products #(r, t) cause variations in the 
perfusion rate ry(r, t) of transport agent at the same point only. Besides, it turns 
out that under these conditions the response of the supplying network keeps the 
concentration of life activity products rigorously inside the vital interval. In other 
words, the field 0(r,i) cannot go beyond the interval (0, 9 C ) for a long time. 

We consider a perturbation 60(r,t) of the field 6(r,t) whose characteristic spatial 
scale Ce meets the inequality 
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(4.1) 



C e > I 



so we may ignore nonuniformities of the field 9(r,t) as well as the field rj{r,t) on 
spatial scales about the elementary domain size In- 

4.1. Linear model for branch response. Perfect self-regulation. Let us 

find the particular form of relationship (2.14) for the function 4>(9) of the form 



(4.2) 




o < e < 

6>6 C 



This dependence <fi ( 9) as well as the corresponding behavior of the f uncti on f (x) is 



displayed in Fig. 2.2 by dotted lines. In particular, as follows from ( 2.13 ) f ld (x) = 1 
for < x < 1 and f ld (x) is undetermined at x = 0. 

Besides, we adopt the following additional assumptions. First, the delay time t„ 
of the branch response is set equal for all the branches, i.e. r„ = r. Second, we will 
ignore the term r^-lnJj in equation ( |3.10| ). This term actually describes dynamic 
interaction between additional potential sources and has no substantial effect at least 
in the following two limit cases. When r < re, where Tg is the characteristic time of 
variations in the living medium state, the network response will be quasistationary, 
and the term In Ji along with the first one in ( 3.10 ) can be ignored. When r > 
and, in addition, the generation rate q is not too high the value of £iT-^ In Ji should 
be negligibly small in comparison with lnej. Indeed, in this case variations of the 
field 9(r, t) are too fast for the supplying network response to follow them and the 
values of {£j} should be small. Since time variations of the transport agent flow Jj 
are due to the action of the additional potential sources {£j} we may write Ji oc {£j} 
and, thus, the term e^T^lnJ^ is of second order in {ef\. Besides, the inequality 
£i < JiRn is considered to be fulfilled in advance (below in the present section we 
will justify this assumption) and so / ld (l — Ej/jJiB^. )) = 1. 
Under these conditions we can rewrite equation ( |3 . 1 ) as 



(4.3) 



where we have also taken into account relation ( |3.12 ). Then acting by the operator 



(1+ tJj) on equation ( ft.9| ) and using (4.3) we reduce it to the following 



(4.4) r | + , 



Vo 



T I dv 'Y.W E e^X^r'Mr'Mr'). 



Here rjo is the perfusion rate of transport agent through the living medium containing 
no life activity products, 9(r,t) = 0. It should be pointed out that the value 770 may 
be regarded as a certain phcnomenological parameter in context of the self-regulation 
problem. 



By virtue of inequality (|4.l| ) we can simplify equation (L4) in the following way. 
In the case under consideration there is a spatial scale I* such that (jv C f < £e on 
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which the fields rj(r,t), 0(r,t) are practically constant. Let us choose a branch level 
n* 3> 1 for which l n » w £* and N — n* 3> 1. At the given level it is possible not to 
distinguish between the true perfusion rate T)(r,t) and the perfusion rate T)*(r, t) ~ 
r/(r, i) averaged over the size l n * of the fundamental domains M„* . So we may specify 
this averaging, for example, in a manner as it is done in Appendix JbL namely, 77* (r, t) = 
Pn* {v( r i t)}- Th en a cting by the operator P n * on equation (|4.4|) and taking into 
account identity ( B.3 ) we convert it to the same form within the replacement 



w E - E e *( r ) A 



The latter, however, enables us to make use of identity (BA) and to reduce (4.4) to 



(4.5) 



dr] 



T f dv 'w E Qi{r)Qi{r'WW). 



Whence it immediately follows that 



(4.6) 



'TO 



because the fields 6(r, t), f)(r,t) are practically constant on scales of order l n *. 

Equation (L6) is the desired relationship between the perfusion rate rj(r,t) of 
transport agent and the concentration 8(r, t) of life activity products. 

4.2. Characteristics of perfect self-regulation. The local form of equa- 



tion (4.6) demonstrates that under the given conditions the perfusion rate r](r,t) 
is determined by variations in the field 9(r, t) at the same point only. Thus, the co- 
operative response of all the branches is so self-consistent that the supplying network 
delivers "nutrients" only to the living medium points that "ask" for this. In this case 
no unexpected changes in the "nutrient" delivery interfering with the living medium 
activities occur, i.e. the supplying network functions perfectly. 

Another characteristics of perfect self-regulation is the fact that the field 8(r, t) 
cannot go beyond the vital interval [0, 9 C ] for a long time. This is due to an unbounded 
increase in the perfusion rate i](r,t) that occurs when the field 6(r,i) locally exceeds 
the value 9 C (see equation ( }4.6| )). Typical dynamics of the fields 0(r,t), rj(r,t) demon- 
strating such behavior is shown in Fig. 4.1 for the one-dimensional living medium. 

Concluding this subsection we remind that when deriving equation ( J4.6| ) we have 
used the inequality e,; < JiR„ without justification. Its validity can be shown in the 



following way. Let us consider the potential drop Ai^ = P? v 
i whi ch according to (3.7) is equal to AP; = JiR^ — £j. Then taking into account 



( p. Ill) , Q4J), and 0461) we get 



P* ut across a branch 



(4.7) 



dAP, 
dt 



AR = R° 



1 In l dr9j(r). 
on 



The right-hand side of this equation is a constant independent of variations in the 
living medium state. At the initial time the field 9(r,t) = and so APi > 0, thus 
APi should be positive at any time moment later. 
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Fig. 4.1. Typical dynamics of the fields 8(r,t), r](r,t) under local stepwise variations in the life 
activities of the living medium (q(r,t) is the generation rate of life activity products). 



Besides, equation (4.7) demonstrates one more characteristic property of the per- 
fect self-regulation. Under the given conditions the potential distribution over the 
supplying network does not depend on variations in the living medium state. Exactly 
this independence ensures that the transport agent flow Ji through a last level branch 
i € Ljy will remain unchanged if the state of the living medium inside the correspond- 
ing elementary domain 971,; 6 docs not alter. In fact, for this branch the value 
Oi = 0(r,t) (r 6 so, the "resistivity" Ri(6i) of branch i will be also constant 

under such conditions. Consequently, the transport agent flow Ji — APi/Ri(9i) as 
well as the perfusion rate r](r,t) = Ji/l^ in this domain will be independent of the 
living medium state at other points. In other words, it is the constancy of the poten- 
tial distribution that ensures the perfusion rate ry(r, t) being controlled by the living 
medium state at the same point. 

5. Model generalizations. 

5.1. Dichotomic supplying network. In the previous sections we have con- 
sidered the supplying network whose architectonics and embedding into the living 
medium are specified by g-fold spliting of branches with g = 2 d , where d is the living 
medium dimension. However, looking attentively through the previous analysis we see 
that it has been used, in fact, only the following properties of the supplying network 
architectonics. Each branch i supplies (or drains) as a whole a certain domain 9JT 
and all the branches of smaller length and connected with the branch i supply (or 
drain) subdomains of the domain 371*. The whole collection {97T}i e L„ of the domains 
corresponding to branches {i S L„} of the same level L„ compose together the living 



medim 9Jt and are mutually disjoint (see (2.1)). Besides, if {i S L^} is the whole 
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Fig. 5.1. Fragments of the supplying network architectonics for the two- and three-dimensional 
living medium. 



collection of branches of level n that are formed by splitting of the branch j (of level 
m < n) then all the fundamental domains {TV} together make up the domain SOP 
(see (f?~^)). None of the other properties has been used. 

In particular, the factor g is not resticted by the value of 2 d . Therefore the ob- 
tained results hold for any other implementation of the architectonics of the supplying 



network and its embedding into the living medium that meets conditioins (2.1), (|2.2|) 



For example, dichotomic supplying networks demonstrated in Fig. 5.1 obey these 
conditions. 

5.2. Irregular supplying network. In the previous sections we have analyzed 
the self-regulation problem assuming the supplying network regular in architectonics 
and involving a large number of levelvs. So it could give the impression that these 
assumptions are the necessary conditions for self-regulation to be perfect. Therefore 
in the present section we demonstrate that the perfect self-regulation can occur in sys- 
tems with supplying networks of arbitrary architectonics, at least, when the transient 
processes are ignorable. 

As before we consider a supplying network consisting of supplying and draining 
beds which are mirror image of each other and connected through terminal branches 



(Fig. 5.2). The transport agent does not interact with the living medium until it 
reaches the terminal branches where the transport agent and the living medium ex- 
change "nutrients" and life activity products. Then during motion through the drain- 
ing bed the transport agent does not interact with the living medium again. Transport 
agent flow is assumed to be governed by the same regu larities as stated in §||. 

So for each node b of the draining bed (Fig. |5.3| a) we can write the following 
conservation laws 



(5.i) ^2 j i = J2 Ji and Y, 9 ^ = ( out)b^^, 

z£b i£b iEb iEb 

where the sums X^eb an< ^ Sieb run over au " the branches i £ b leading to or from 
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terminal branches 




Fig. 5.3. Elements of the irregular supplying network. 



this node, Ji and 6i are the transport agent flow and the concentration of life activity 
products in the branche i going in or out of the node b. Besides, for all the branches 
going out of one node b the concentration (0 O ut)& of life activity products is assumed 
to be the same. The branch orientation is chosen according to the direction of the 
transport agent motion. The potential distribution {Pb} over the draining bed nodes 
and the patten {Ji} of transport agent flow are related by the expressions 



(5-2) J i R i (6 i )=I&-I* at 'gAP i , 

where P? n and P* ut are the potentials ascribed to the terminal nodes b\ n and b l out of 



the branch i as shown in Fig. 5.3b. The dependence of the "resistance" Ri{9i) of the 
branch i on the value 9i describes the living system response to variations in the living 
medium state. In the present section we analyze only the case corresponding to the 
perfect self-regulation which is characterized by the relationship (cf. ([2.12])) 



(5.3) Riifii) 



Uc 

if 6i > 1 



COOPERATIVE MECHANISM OF SELF-REGULATION 



19 



where R® is the "resistance" of the branch i at 9f = 0. For the terminal branches 
{i } term of the draining bed the quantities {#i} tcrm are directly specified by the living 
medim state and, so, may be treated as predetermined when analyzing the distribution 
of transport agent flow over the draining bed. Besides, due to the mirror symmetry of 
the supplying and draining beds, we can regard the potential P° ut at the exit node of 
the draining bed (the latter node of the draining bed stem) and the potential P^ Im 
at all the entrance nodes (the former nodes of the terminal branches) as fixed values. 

In order to show that under such condition the supplying network response is 
perfect it is sufficient to demonstrate that the potential distribution over the draining 
bed nodes is constant, i.e. does not depends on variations in the field 9(r,t). Indeed, 
in this case the perfusion rate ry(r, t) of transport agent inside a domain 97t* orm drained 
as a whole by the terminal branch i term is 



~ J i° rm AP* ™ const 

where y. term is the volume of the domain 97t* crm . Whence it follows that the perfu- 
sion rate ?y(r, t) is practically controlled by the field 9(v, t) taken at the same point 
r S97l* crm because the value of 9\ crm is determined directly by the value of the field 
0(r,i) in the domain 9Jtf rm . 

Let us fix a potential distributi on { Ps } th at occurs at initial time. In general 



the system of Kirchhoff's equations (5.1), (5.2) should have a single solution. So to 
prove the independece of the potential distribution {Pf,} from variations in the field 
0(r,t) we may to show that this system possesses a solution provided we have fixed 
{P b =P b }- 



Let us consider the system of equations (5.1), (5.2) for an arbitrary node b 
(Fig. |5.3| a) which is formed by g^ n branches going into it and g^ ut branches leav- 
ing this node. At first, we will treat the values {0i} m ascribed to the branches going 
into the node b as given quantities. Then the variables to be found by solving this 
system are as follows: { Ji} ln , {Ji} ont , and (# ou t)b (the subsripts in and out mean that 
the corresponding quatities relate to the branches going into or out of the node b). 

The total number of these quantities is equal to gf n + g^ ut + 1, whereas the total 
number of the given equations is g^ n + g% ut + 2. Therefore this system seems to be 
overdetermined. In fact, however, equations ( |5.1| ), ( |5.2| ) are linearly dependent. 

To justify this statement, we, first, note that under the steady-state conditions 



none of the quantities {9i} can exceed the value 9 C . Indeed, according to (5.1), for 
any node b the value of (0 O ut)b < nrax{#i}i n . So for a branch i the value of 9i might 
exceed 9 C only if there were a path Pj on the draining bed leading from the branch 
i to one of the terminal branches such that for all braches {j} belonging to it the 
values 9j > 9 C . In this follows from ([Tl) and (|5.3D, the total "resistance" of 



the path P, were equal to zero. So the potential P* nt at the node bj into which the 
branch i goes would coinside with the potential at the draining bed entrance, P* ut = 
ptcmi j n addition the total transport agent flow J2iTb Ji g° m g through the node bi 
would be equal to the transport agent flow along the path Pj except for cases when 
there are another similar paths P^. This in turn would cause the value of (9 out )b i to 
become greater than (or equal to) 9 C and so on. In this way we would reach the stem 
of the draining bed and obtain that the potential P^ ut at the exit of the draining bed 
has to be equal to the portential P^ lm at its entrance. Thus all the values 9i do not 
exceed 9 C . 
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Then let us sum equations (5.2) with the wight 1/R® over all branches relating 



to the node b. In this way taking into account (5.1) and (5.2) we get 

111 pi p out p pi 

(5.4) E^V^ = E b "'" 



Ri ft R i 



Since the initial portential distribution {Pf,} has been established by the existing 



transport agent flow this equality is fulfilled. So equations (|5Tl|), (5.2) are linearly 
dependent and can be reduced to a system of (gf n + g% ut + 1) equations which is not 
overdetermined and possesses a solution. 

Above we have regarded the quantities {9i} as given beforehand. If a node b is one 
of the latter nodes of the terminal branches this is justified because the values {9i} at 
such branches are di rect ly det ermined by the living medium state. Then, sovling the 



system of equations (5.1)-( |5.2| ) for these nodes we find the values {(# ut)b} which form 
the collection of quatities {0i}i nt for nodes of lower levels. Repeating this proceedures 
we get the draining bed stem. 

I n th is way we get the conclusion that there is a solution of the full set of equa- 
tion (|5.l| )-(j5~i^) for a fixed potential distribution {Pb} when the state of the living 
medium varies in time. The latter demonstrates that under the adopted assumptions 
the response of the supplying network is perfect. 



6. Essence of perfect self-regulation. The fact that equation (4.6) governing 
the dynamics of transport agent perfusion can be of such a simple form is a surprise 
because it does not directly contain any information on the complex architectonics of 
the supplying network. This fact is actually one of the main results obtained in the 
present paper. So we, now, discuss the qualitative physical essence of haw the perfect 
regulation can occur. 

Let us assume that in a small domain (Fig. |6.l| ) the field 9(r, t) increases 
and comes close to the boundary C of the vital interval. In order to smother this 
dangerous growth the system should increase the perfusion rate r/(r,t) in the domain 
0. For this purpose all the "resistances" {P,}p along the path P connecting the 
stems of the supplying and draining beds and going through the domain have to 
decrease substantially. As a result, the transport agent flow along the path P and, 
consequently, the perfusion rate in the domain will grow. Higher values of the 
variables {#i}p along the path P (on the draining bed) bear the information required 
of such a response. 

Variations in the field 0(r,t) located in the domain can, in principle, give rise 
to alterations of the perfusion rate at points external to the domain because of 
variations in the transport agent flow through large branches. However there is a 
difference between the external and internal points of the domain from the stand- 
point of this transport agent flow redistribution. For an external point r ^ on the 
path P there are two nodes b£ an d belonging to the supplying and draining beds, 
respectively, at which a similar path P r leading to the given point r, at first, diverges 



from the path P and, then, they converge (Fig. 5.1). Due to the mirror symmetry of 
the two beds we can confine our consideration to the supplying bed only. The node 
b£ divides the path P on this bed into two parts. One of them leads from the stem to 
the node 6* an d coincides with the path P r . The other goes from the node b* to the 
domain and does not belong to the path P r . Therefore, on one hand, decrease in the 
"resistances" of the branches forming the former part tends to increase the transport 
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Fig. 6.1. Illustration of the mechanism by which the nonlocal dependence of the perfusion rate 
n{9} on the field 0(r) at distant points is suppressed. 



agent flow along the path P and, thus, through small branches of the path P r that 
directly control the perfusion rate in the vicinity of the point r. On the other hand, 
the less are the "resistances" of the branches making up the latter part of the path 
IP, the larger is the fraction of transport agent flow that is directed along the path 
P at the node These effects are opposite in sign with respect to increase in the 



perfusion rate at the point r as illustrated in Fig. 6.1 and can compensate each other 



In the present paper we have actually found the conditions of such compensation. 



7. Conclusion. In the present paper we have shown that there is a specific co- 
operative mechanism of self-regulation by which a living system can respond perfectly 
to changes in the environment, i.e. react without interference with itself. This mecha- 
nism is based on the corresponding self-processing of information and the cooperative 
effect caused by individual response of the supplying network elements. The existence 
of large hierarchical systems in nature and their capability for adapting to changes in 
the environment points to the fact that they should function perfectly, at least, at first 
approximation. For example, in living tissue in order to reduce effects of nonideality 
arterial and venous beds contain a system of anastomoses, i.e. vessels joining arteries 
or veins of the same level jlj. In other words, the vascular network is organized in 
such a manner that its function be as perfect as possible. 

So we think the obtained results will be useful in analysis of real living systems and 
can form the basis of the following mathematical modelling of complex self-regulation 
processes. 
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Appendix A. The Green matrix Ay. 

By definition, the Green matrix Ay gives the flow pattern { J; = Ay : j is fixed} on 
the homogeneous draining bed (Fig. A.l) provided the branch j is fixed and contains 
a single potential source of pow er Ej = 1. (For simplicity of dra wing we de pict the 



binary draining bed in Fig. AT.) In other words, as results from (|2.5|) and (g^D the 
Green matrix Ay is the solution of the following system of the conservation equations 
at the nodes {b} 




Fig. A.l. Homogeneous draining bed with a single potential source of unit power, Sj = 1, 
located at the branch j . (For simplicity is drawn the binary draining bed, g = 2.) 



(A.l) £A^ = (A^)L, 

ieb 

and the potential equations 



(A.2) 

Here the sum runs over all the branches ieb going into the node b and 6ij is the 
Kronecker symbol (Sij = 1 if the branches i, j coincide with each other and 5ij = 0, 
otherwise) . 

In order to specify the Green matrix we classify all the possible pairs of 
branches into two groups. The first group comprises all pairs {i,j}+ that can 

be joined by a path on the draining bed directed either from higher to lower levels or 



vice versa, i.e. by a path of constant direction. Such a path is depicted in Fig. A.l 
by the dotted line P + . The second group involves the pairs that can b e join ed 

by a path whose direction changes at a certain node by as shown in Fig. A.l by 



the dotted curve P . In other words, this path initially goes, for example, from the 
branch j towards the stem until it reaches the node b,j (in Fig. A.l it is the node 



bji- ) and then goes towards the branch i in the opposite direction. To the node by 
we also ascribe the level number of branches going into it. The given classification 
allows us to state the following. 

Proposition A.l. Let us introduce the quatities {Z(n)} (n — 0, 1, . . . , N) such 
that: 
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N 



(A.3) Z(n) = £ p(m), 

m—n 

where the function pirn) is defined by formula (2.11). Then at leading order in the 



small parameter p(n)jZ(n) the Green matrix Aij is specified by the expression for the 
Pairs 

(A.4a) A i:j = Tr^W - !" if Ui ^ n i> 



and for the {i,j}~ pairs 



Ro Z(rii) 

(A.4b) A --i ^zW } if U ^ n - 



(A.5) Aij = - - g ™-n*-™ p{nii) 



(.g-l)i? y [Z( nij )r 

where n t is the level number of branch i. 

The given proposition is the main result of the present Appendix and the remain- 
ing part will be devoted to its substantiation. 

Proof. Let us, first, prove this statement in the case when the branch j is the 
stem io- The values of Au are the same for all the branches {i : i € L ni } belonging 



to one level L rai . So from (A.l) we get 



(A.6) 



A, 



'A,- 



because level L„ involves exactly g n branches. Then choosing any path P on the 
draining bed leading form the stem i to a branch ijy of the last level and summing 
equations (A. 2) along this path we obtain 



(A.7) 



} AjjR^. = A loio R o Z(0) — 1, 



where we have taken into account expressions (2.11) and (A.6). Whence we immedi- 
ately find the expression for Ai a i a and from (AJ3) get the desired formula 



(A.8) 



A., 



1 



RaZ{0) 



which gives the same result as does formulae (A.4) (namely, (A.4b) for 



because all the pairs {iio} belong to the first group {ij}+. 



0) 



To find Aij for a branch j whose level number n,- > 1 we consider a path on the 



network shown in Fig. A.l by the dotted line that goes from the given branch j to the 
stem ig. This path divides the draining bed into the path P° itself and disjoint sub- 
trees whose stems are connected with the path P° through its nodes. The values of A^- 
for all branches of one level belonging to the same subtree are equal. This allows us to 
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Fig. A. 2. T7ie bloc k repr esentation (upper graph) of the homogeneous draining bed equivalent to 



the graph shown in Fig. A.l and the block structure (lower graph). (A block unites the draining bed 
subtrees characterized by identical transport agent flow distribution and r^ T is the total "resistance" 
of the corresponding subtree.) 



transform the graph shown in Fig. A.l into one shown in Fig. A.2 , where the path Fj 
is represented by the sequence of branches designated by {j, j n ._i,j n ._2, . . . io}- 
Here the symbol j n stands for the branch of level n belonging to this path, with j n . and 
jo indicating the same as j and io, respectively. The kinetic coefficients corresponding 
to these branches are {R n ., R n ._i, R n -2> •••> Ri> Rq}- These branches at the ter- 
minal points (nodes bound up with their entrances and exits) {b* , b nj , b nj _i, . . . , bi} 
are connected with the blocks {B* .,M nj ,M nj -i, . . .,Bi} of identical subtrees. These 
blocks except for the former one B*^. (connected with the node b* ) involve (g — 1) 
subtrees, the block B*. contains g subtrees. The transport agent flow spreads over 
such subtrees uniformly. Therefore each block, for example, B n connected with a 
node b n can be treated as a single element characterized by the kinetic coefficient 
( "resistance" ) 



(A.9) 



In 



Here /„ is the total flow of transport agent through the given blook, P$ n is the 
potential at the node b n induced by the potential source ej = 1, and for the sake of 
simplicity we have set the potential at the entrances of the la st le vel branches equal 
to zero, P£ = 0. As is seen from the diagram shown in Fig. A. 2 (lower graph), for 



1 < n < rij the value of r, 



is specified by the expression (see also Fig. |A.2|) 
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1 



1 



-iL, 



1 



(.9-1) " 9-1 " (.9-1)5 

N 



1 

W^)9 



R n 4-1 + 



2 ±L n+2 



(A.10) 



1 Y,9- (p - n) K 



(.9 - 1) ^ 



because level p > n contains (g — l)g p n identical branches belonging to this block. 
Expressions (2.11) and (A.3) enable us to rewrite (A. 10) as 



(A.11) 



1 



(.9 - 1) 



g n R Z(n). 



In a similar way we obtain the expression for the kinetic coefficient of the last block 
connected with the path P° through the node b* 



(A.12) r nj+1 =g n 'R Z(n j + l). 

Let us introduce the quantities 

J n = A jnj and I n = }^ A 3^j 

being the transport agent flow through a branch j n and the total transport agent 
flow in the subtree block B„, respectively, which are induced by the potential source 
Ej = 1. In the last expression the sum runs over all the stems of the subtrees forming 
the block M n . Then the given reduction of the initial graph to the subtree blocks 
(Fig. A. 2) enables us to rewrite the Kirchhoff equations (A.l), (A. 2) in the form for 
1 < n < n,: 



(A.13) 

and for 2 < n < n 7 - 



J\ n Jn — 1 



(A.14) 



In r n — Jn—lR n — l ^n— iTn— 1- 



These equations are completed by the following expressions bound up with the first 
and the last elements of the given network: 



(A.15) 



hn = JqRq 



(A.16) 



J nj [r nj+1 + R n] ]+I nj r nj =1 



Due to p(n) being a smooth function of n the ratio p(n)/Z(n) may be treated as 
a small va lue fo r N — n 3> 1. The latter allows us to find directly the solution of 
equations ( A.13 )-( A.16 ), which is matter of the next step. 
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Co ntinu ous s olutio n of equations ( A. 13 )— ( A..16| ) Takin g into account ( [2.11 ), 
(A. 3), ( A. 11 ), and ( A. 12 ) we can rewrite equation ( A.14 )-( A. 16] ) in terms of 



(A.17) glnZ(n) = {g- l)J n -ip(n - 1) + I n -iZ(n - 1) 

for 2 < n < rij and 



(A.18) 



gI x Z(l) = (g- 1) J , 



(A.19) 



-I, 



(.9-1) nj R Z( nj ) 



9 



where, in addition, we have used the identity Z(n + 1) + p(n) = Z{n) and the equality 
p(0) = 1. 



Let us analyze, at first, the solution of equations ( A. 13 ) and ( A.17 ) when n ~ 1. 
For such values of n the function p(n) as well as Z(n) may be considered as being 
constant, p(n) = p(0) = 1 and Z(n) = Z(Q) 1 due to smoothness of the function 
p(n). We seek the solution of these equations in the form 



(A.20) J n =A+Cf+A_C, 

where are some constants. Substituting (A.20) into ( A.13| ) we get 



(A.21) 



in = A + d n - 1} (c + i) + A_d n - 1} (c- - 1) 



Then the substitution of this expression into (A.17) shows us that the constants £+j C- 
are the roots of the equation 



(A.22) 



(C-i)(,C-D = |^C. 



Whence it follows that to the first order in the small parameter 1/Z(0): 



(A.23) 



1 



Z(0) 



, C- 



g 



l 



Z(0) 



To the same order in l/Z(0) from (A. 13) and ( A.lq ) we find that the constants A +1 A- 
are related by the expression 



(A.24) 



A_ _ 1 
17 ~ (.9 - WO) 



< 1. 



According to ( A.23 ) the first terms on the right-hand side of ( A.2C ) and ( A.21 ) are 
increasing fu nction s of n, whereas the second ones are decreasing functions. Therefore, 
by virtue of ( A.24 ) in the two expressions the former term is substantially greater than 
the latter for all n *~ 1. So, we may ignore the second term and regard the first one 
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as a function of the continuous variable n because of (C+ - 1) < 1. This allows us to 
set 



(A.25) 



'J 7i Jn—1~ 



dJn 

dn 



in equation ( A.13| ) and I n -\ ~ I n m ( A..17 ), and also to rewrite ( A.2C ) in terms of 



(A.26) J„~A + exp{^}. 

In this case the system of equations ( A.13| ) and JA.17 ) is reduced to the equation 



(A.27) 



dJn 

dn 



Z(0) 



and function ( A.27 ) is its general solution. So, the "b ounda ry cond ition" ( A. IS ) is 
responsible only for the existence of the second terms in ( A. 20 ), ( A. 21 ) and, therefore, 
can be ignored. 

As follows from equations ( |A.13| ), ( A. 14 ) (or ( A. 17 )) their general solution con- 
tains two arbitrary constants, for example, Jo, J±. Indeed, all the other values of 
I n and J n can b e foun d by ite ration. The same situation is met with respect to the 
general sol ution ( A.2C ), ( A. 21 ) of this e quatio ns for n ~ 1. So there is no solution of 
equations ( A. 13 ), ( A. 14 ) different from (|A.2C), (A. 21 ) for such values of n. Th ereby , 
the i nfluen ce of the "boundary condition" ( A. 15 ) on the solution of equations ( A. 13 ) 
and (A. 14) is ignorable and we may seek this solution in the class of smooth functions 
of the continuous variable n. In this case it contains only one arbitrary constant 
specified by the "boundary condition" (A. 16) (or (A. 19)). 

Within the framework of the adopted assumptions the quantities J n and I n as 
well as p(n) and Z(n) are related by the expressions 



(A.28) 



dJn 

dn ' 



(A.29) 



dZ(n) 
dn 



and the system of equations ( A.13j ), ( |A.17 ), ( A.19| ) is reduces to the following equation 



(A.30) 

subject to the boundary condition 



dJn y,., T dZ(n) 



(AM) 



Jn, = 9~ 



RoZ(nj 



In obtaining these expressions we have also taken into account relation ( A. 13 ) and 
the inequalities \I n +i — In\ "C I n and \Z(n + 1) — Z(n)\ <C Z(n). Besides, we have 
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ignored the second term on the left-h and side of ( A. 19 ) because of /„ . = J, u ~ i — J n <C 
J n . . The solution of equation ( A. 30 ) meeting the boundary condition ( A. 31 ) is of the 
form 



(A.32) 



1 

RaZ(n) 



and, thus, from ( A.2§| ) and ( A.29| ) we obtain 



(A.33) 



In = 9' 



p(n) 



For a branch i located on the path P°, i.e. i £ {j, j nj -i, jn. 



,ji,i } the pair 



{ij} belongs to the first group and Ay = J ni (rii is the level number of the branch i). 



Whence we directly get formula (A.4a). If a branch i belongs to the last block 



of subtrees it level number n.; < n.j and the pair {ij} also belongs to the first group. 
Since for a fixed value of the total number of such branches is g ni ~ n i and the 
transport agent flow (Ay ) is the same for all of them, we ca n write Ay = J n J g ni _rlj 1 , 
which together with expression ( A.32 ) gives equality ( A. 4b ). 

Now let us consider a branch i located in one of the blocks {B„ , B„ i, . . . , Bi}, 



e.g. in a block 



For this branch the pair {ij} belongs to the second group and 



the node b at which the block B n „ joins the path P° plays the role of the node by 
of the pair {ij}-. For a fixed level rn > nj the total number of such branches is 
(9 ~ l)g ni ~ m and the total transport agent flow going through this bl ock is I nb . So 
in this c ase A y = —I nB /[{g — whence taking into account ( A.33 ) we get 

formula ( |A.5| ). In this expression the sign "— " means that in the given branch i the 
transport agent flow Ay caused by the additional potential source £j — 1 is directed 
from lower to higher levels. In other words, this component of the transport agent flow 
through the branch i and the total transport agent flow J,; induced by the collective 
action of all the potential sources have opposite directions. 

In this way we have considered all the branches, so, proven the stated proposition. 

□ 

Appendix B. The embedding identities. 

Let us specify the averaging procedure V n ■ \P(r) — > W(r) of a function ^(r) over 
spatial scales of the fundamental domains {9Jt n } of level n by the expression 



(B.l) P n {y(r)} d ^ J dr'F„(r,r')*(r'), 

wi 

where the kernel 



(B.2) F„(r,r') = - r ^e 4 (r)a i (r') 

and the sum runs over all the branches {i} forming the level L n , for example, of the 
draining bed. Then we can state the followig. 

Proposition B.l. Let L„ anc?L m be the branch collections forming the levels n 
and m, respectively, and n > m. Then 
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(B.3) 



-1^6,(^=^(1^6^,1. 

m ieL m I « ieh n ) 



Proof. We prove identity (B.3) by the direct substitution of (B.2) into the right- 
hand side of (B.3). In this way we get 



m iGL,„ t'elLj, 

Here we have taken into account that, first, 



1 



7 / dr , e i (r')9 i /(r') 



1, if i' e L* 
0, if i'rfLi 



an 



where is the full collection {i'} of branches of level L n belonging to the subtree 
whose stem is the branch i and, thus, £Dtj/ C SEJtj, and, second, 



i'ei>* 



because of the conservation laws (A.l). □ 

In the present paper we consider the model for the supplying network which 
envolves large but fininte number of hierarchy levels. In other words, we regard the 
total number N of levels as a fixed (but large) value. In this case we can state the 
following. 

Proposition B.2. At leading order in the small parameter p{ni)/Z(yii) for al- 
most all points r S 9JI it is fulfilled the identity 



(B.4) A 4 X j e i (r)=e i (r), 

provided the level number Hi of the branch i meets the inequality 1 -C Ui <C N. 



In this context the term "leading order" means, in addition, that in (B.4) we 
ignore any summand caused by variations of the quantities p(n) and Z(n) when the 
argument n changes on scale of order unity. This is due to such summands cannot be 
analyzed in the framework of the adopted continuous approximation of Kirchhoff 's 
equations. 

Let us prove this assertion. 

Proof. For a finite number of hierarchy levels the volumentric measure of the 
interface made up of the boundaries of all the fundamental domains is equal to zero. 



So in proving ( B.4 ) we may ignore all the points belonging to this interface. Therefore 
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OTjv-i(r) m N (r) m N . 2 {r) 




b* 

Fig. B.l. Illustration of constructing the domain collection {9Jln(r)} and the collection of the 
corresponding branches {i n (r)} forming the path P(r) (solid line) on the draining bed (dotted lines). 
(For simplicity is shown the two-dimensional living medium.) 



let r be a point of the living medium domain, r G 9Jt, that does not belong to the 
boundary of any fundamental domain. For each level number n this point r is located 
inside just one fundamental domain of level n which will be labelled as 07l n (r). In this 
way we can build up the collection {iUln (r)| (n = 0, 1, . . . , N) of all the fundamental 



domains containing the point r (Fig. B.l ). These domains in turn specify the sequence 
{i n (r)} — {i , ii(r), ia(r), . . • , iAr(r)} of connected branches whose n-th term is the 
branch i n (r) of level n contained in the domain 97l„(r). The branches {i n (r)} form 
a continuous path P(r) on the draining bed 9T d that leads from the stem i to the 
branch lAr(r) of the last level. 

In formula (B.4) the functions {8j(r)} differ from zero for the branches {i n (r)} 
only, thus 



(B.5) 



r . . dcf ^— \ 

A 4 (r) = 



^ A,, /<■'' H,fr; 



E 

j£P(r; 



A R° 



Let us, first, consider a branch i being one of the branches {i n (r)} whose level num- 
ber rii meets the inequality N — rii 3> 1. For the given branch all the pairs {i, i n (j)} 
belong to the first group specified in Appendi x |A . Then taking into account expres- 
sions (2.11), (|A.3|), and (A. 4) we can rewrite (B.5) as 



A,(r) = J2 A «< + E A «< 

jEP(r) jEP(r) 

1 n m-i , , 

~E PK-) + E^ (n< -" 3) |r4 

Z(rii) ' '-^ Zirij) 

<■ l > n j= ni n j= v ■* > 
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and, thus, 
(B.6) 



A i (r) = l + £> 
P =i 



at leading order in the small parameter p{rii)/Z(rii) for N — rii 3> 1. 

Let us, now, consider a branch i not belonging to the branch collection {i„(r)} 
(in Fig. B.l it is the branch i*). In this case there are a single node b (node b* in 
Fig. BT) on the path P(r) and a single subtree that contains the branch i and whose 
stem joins the path path P(r) through the node b. The node b devides the path P(r) 
into two parts, one going from the last level branch ijv(r) to the node b and the other 
leading from the node b to the stem of the draining bed. Let us ascribe to the node 
b the level number rib of branches going into it. Then we can say that the former 
part of the path P(r) involves the branches {i„(r)}_ whose level numbers n > rib- 
For such branches the pairs {«,i„(r)}_ belong to the second group (see Appendix |X|) 
and rib — n u n (r)- The latter part is made up of the branches {i n (r)} + whose level 
numbers n < rib and f° r these branches the pairs {i, i n (r)} + belong to the first group. 
In the given case taking into account expressions ( |2.11 ), (A. 3), ( |A.4| ), and (A. 5) we 
can rewrite (B.5) as 



j'eP(r) jeP(r) 



1 



(.9-1) 



9 



(.9-1)- 



[Z(nbW 

P( n b) , 
Z(n b ) 



E P(*i) + E 5" 

nj>n b rij<n b 
rib 



P =i 



Z(nj) 



g^-my g -P 



P(n b -P) 
Z{n b -p) 



Whence treating the small ratio p(rib)/Z(rib) as a smooth function of the continuous 
variable n we get for rib 3> 1 



(B.7) 



A*(r) 



d 

dn 



pin) 



Z(n) 



Era" 



at leading order in this small parameter. Therefore we may state that Aj(r) = o( §fe4 ) 
because in the given case rib < n i an d so either rib ~ i*i or n, - n t > 1 and for such 
values of rib the magnitude of Aj(r) is of exponential smallness: Aj(r) ~ xhe 
given nonzero value of Aj(r) is actually due to variations of the function p(n)/Z(n) 
as the argument n changes on scales of order unity, so in this case we have to set 
Mr) = 0. 

Summarizing the aforesaid we see that if the point r is located inside the fun- 
damental domain dJli corresponding to the branch i, i.e. Sj(r) = 1 then the branch 
i belongs to the path P(r) and the value Aj(r) = 1 at leading order in p(n)/Z(n). 
Otherwise, when r ^ dJli and so i $ P(r) the value Oj(r) = and we have to set 
Aj(r) = 0. This substantiates identity (BA). □ 
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